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A PRACTICAL METHOD OF DETERMINING WATER CURRENT 
VELOCITIES AND DIFFUSION COEFFICIENTS IN 
COASTAL WATERS BY REMOTE SENSING TECHNIQUES 

by 

Wesley P. James 


INTRODUCTION 


Population and industrial centers along the 
coastal areas are generating large volumes of waste 
materials which present a threat to the environmental 
quality of the area. The environmental characteristics 
of the coastal area can be expected to generate a large 
population growth in the near future which will compound 
the pollution problems of the area. 

Restricted mixing and circulation character- 
istics of the estuaries and nearshore waters limit the 
natural capacity to handle concentrated waste loadings 
from population centers. The major factors that deter- 
mine the mixing characteristics in the coastal waters 
are (1) wind, (2) tidal range, (3) topography, f4) fresh 
water inflow and (5) density gradients. The salinity or 
density gradients are functions of the fresh and salt 
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water inflows and the circulation within the estuary. 
Energy must be provided to break down these gradients 
to allow complete mixing. 

Diffusion coefficients and current velocities 
in the receiving waters are essential to waste disposal 
studies. Mathematical models of the coastal area require 
this information to predict the fate of future waste 
loadings. Since the coastal waters cover large areas 
and the diffusion coe f “ ‘ cients ani current velocities 
vary with both space and time, remote sensing provides 
a practical method of gathering this information. 

Field work begins by dropping dye markers 
from an aircraft at selected locations throughout the 
receiving waters. The movement and spread of the re- 
sulting dye patches are measured remotely by two flights 
over ;he area. The movement of the patch between flights 
gives the current velocity while the change in size of 
the patch is used to determine diffusion coefficients. 
While the procedures described herein utilize aerial 
photography, the method can be adopted to automatic 
identification of dye patches and automatic processing 
with mult i spectral scanning imagers (Elliason, 1971) . 

By using the ratio of light returned in two bands, one 
band in the region of maximum absorption and the second 
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band iri the region cf maximum re-emission of the dye, 
the dye patch can be identified and computation com- 
pleted. 

DIFFUSION’ STUDIES 

Numerous investigators have employed solu- 
tions to the diffusion equations for the estimation of 
waste concentrations in a waste field. If the scale of 
the current eddies is much smaller than the dimensions 
of the waste field, then the Fickian form of diffusion 
equation can be applied. The basic equation is: 


3W 

2T 


(D 

3Y 1 y 9Y' 


— (D — } + - — (D — ) - [— (V W) 
9X 1 x 3X ; 9Z z 9Z' l 9Y ^ Y 1 


* (V W) + (V W)] + S (1) 

ax x 9z z 

where V is velocity, W is waste concentration, D is 
eda/ diffusivity, X, Y, and Z are space coordinates and 
T represents time. The first three terms on the right 
are the diffusion terms, the next three are convection 
terms and S represents the sources and sinks. 

Solutions to the equations have required vari- 
ous assumptions such as steady state condition, no vertical 
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or longitudinal mixing and unidirectional transport 
velocity in the X direction. With these assumptions, 
the equation becomes: 

V = ~ (D — ) + aW (2) 

x 3X 9Y 1 Y dY J K J 

where a is a first order decay constant and aW repre- 
sents a sink or loss in the system. 

Investigators such as Pearson (1955-1967) 
Brooks (1960) and others have reported solutions to 
the diffusion equation for various conditions. Pearson 
points out that for a source with 1) steady unidirec- 
tional current, 2) uniform mixing of the waste ana 
3) continuous unifora flow from the source, solution 
to the diffusion equation is as follows: 


S 


OK 


2 . 35 (1/1) V X 
1 l 

Q 


(3) 


where S Qm is the minimum dilution along axis of waste 
plume at distance X from a point source; D^, is the 
assumed diffusivity in ftVsec; X is the distance 
from source in feet; V x is the average velocity of 
water mass in ft/sec; Q is the waste discharge, MGD 
and d is the assumed mixing depth in feet. 
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Including the decay function for bacterial 
dieaway or disappearance, and expressing the waste con- 
centration in terms of coliform concentration, the 
above expression becomes: 


0.42S QC 0 

MPN = (4) 

d/DyYT exp (zXV x ) 


where MPN is the most probable number of organisms per 
ml on plume centerline at X; Co is the concentration of 
organisms in waste, MPN/ml ; and a is the bacterial dieaway 
(decay) constant, sec -1 . 

Pearson (19f 7) further points out that: "The 

above equations assume a constant eddy diffusivity; cor- 
respondingly, the value of employed must be represen- 
tative of the overall or average scale of the diffusion 
phenomenon". 

Brooks has reported a solution to the diffusion 
equation with a variable coefficient of diffusivity. It 
’s assumed that the diffusivity coefficient, D^., varies 
as the four-thirds power of the scale of the diffusion 
phenomenon. = al 1 * / 3 , where a is a constant. Brooks' 
equation for a line source is as follows: 
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( 5 ) 
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where: Co **' initial coliform concentration 

C ** maximum coliform concentration at time, t 
m 

t = time of travel X/V x 


IS * 12I) y /V x h 
a = decay constant 

D y = eddy diffusivity at source (X a 0) 
b = initial width of sewage field 
Considering the foregoing solutions to the 
diffusion equation, the four characteristics of the 
receiving waters that have a major effect on waste 
concentration are the following: 

1. V , average current speed 

X . 

2. D , eddy diffusivity 

' \ 

3. d, average mixing depth 

4. a, decay or dieaway constant of pollutant 
The Allen Hancock Foundation conducted an in- 
vestigation on the dilution and dispersion of a waste 
field in the sea (1964J. j p their study Rhodamine B 
dye was introduced as 1) slugs from a point source, 

2) a continuous plume from a point source and 3) a con- 
tinuous plume from a volume source. Dye concentration 
was measured with a Turner fluorometer. The mathematical 
models used for analysis of data were statistical models 
based on Gaussian distribution. The basic three-dimen- 
sional model for the dye slug was 
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W(x, y, z, t) = 


exp - 


M 

tt/Z¥ [a x 2 a y 2 o‘ z 2 ] 1 / 2 


2 a x 2 


25 y z 


2 ° z 2 


( 6 ) 


where W is the average concentration of a point X, Y, Z, 
at time t, M is the amount of dye initially discharged 


from an instantaneous point source and a^ 2 


°y 2 > a z 2 are 


the average values of the variances of the concentration 
dis tribution . 

Several of conclusions of the Allen Hancock 
investigation are listed below. 

1. The rate of vertical diffusion can contribute 
significantly to the overall diffusion process 
at wind speeds greater than eight knots and/or 
low water column stability. 

2. The rate of longitudinal and lateral diffusion 
appeared to be influenced by wind speed but not 
by water column stability. 

3. The "4/3 law" relating the lateral coefficient 
of eddy diffusion as a function of average eddy 
scale did not hold in the particular oceanic 
areas studied. 
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Vertical mixing does occur in the waste field 
as well as horizontal mixing. As indicated by Wiegel 
(1964), vertical mixing is difficult to study in the 
laboratory because of limitations of tank size. In 
these studies 'he wind drags the surface water to the 
down wind end of the tank producing a hydraulic head 
which causes a flow in the opposite direction. 

Laboratory studies have indicated that wind 
drag on the water surface produces very little mixing. 
However, when wind generated waves appear, extremely 
rapid mixing occurs us wind waves are rotational in the 
generating area. Masch (1961) conducted a wave study 
in a wave tank and developed the following relationship 
for the coefficient of eddy diffusivity: 

Dy = 0.007-8 (Vs + Qw) 3 - 2 (7) 

where Vs is the surface current and Qw is the water 
particle orbit speed (Qw = H/T, H = significant wave 
height and T = average wave period) . 
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DYE PATCH STUDIES 

The three pictures shown in figure 1 were 
taken July 7, 1969 using panchromatic film type 8401 
with a Wratten 25A filter. While this is not the best 
film- filter combination for observing the dye patch, 
the change in shape of the dye patch can be seen. The 
dye was dropped at 12:19 and the first photo was taken 
at 12:25 from 3000 ft at which time the size of the dye 
field was 160 ft by 40 ft. The photo in figure lb was 
taken at 13:13 from 4000 ft. After 54 minutes from the 
time the dye was dropped, the dye patch had grown to 
approximately 70 ft wide and an overall curved length 
of 1000 ft. The photo in figure lc was taken an hour 
and a half later at 14:43 from 5000 ft. The dye field 
at this time is 2100 ft long and 1300 ft wide. 

The example given in figure 1 was an extreme 
example of elongation, curvature and striation of a dye 
patch. The wind was 5 to 12 knots with a swell height 
of 4 to 6 ft and a water current velocity of 0.4 ft/sec. 
Most dye patches observed have been elongated in a direc- 
tion nearly parallel to that of the water flow. Striations 
and curvature of the dye patch are common; however, in most 
experiments the overall shape of the patch resembles an 
ellipse. 
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c 

At 14:43 from 500(J ft 


Figure 1. Aerial photographs of a dye patch on July 7, 1969. 



The elongation of the dye patdi in the direction 
of flow suggests that dispersion due to a vertical velocity 
gradient can be an important consideration. The upper 
layers of water are first influenced by a change in wind 
velocity or direction and vertical water current gradients 
arc formed. Wind waves and swell both have specific orien- 
tation and may cause diffusion to occur at different rates 
in a horizontal plane. 

The basic diffusion equation given as equation 1 
is reduced to a two dimensional model. 


oW 

3 1 
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3 2 W 
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aW 


( 8 ) 


A solution is 


W(x, y, t) 


W 

max 


exp 



( 9 ) 


where the coordinate axis is assumed to move with the 
dye patch. in the equation W represents the dye concen- 
tration, X and Y are the coordinates from the centroid 
of the dye patch parallel and transverse to the direction 
of flow, D and D are the longitudinal and lateral dif.- 
fusion coefficients, a is a first order decay coefficient 
(which includes the loss to the lower layers due to verti- 
cal diffusion) and o x 2 and a y 2 represent the variances in 
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the X and Y directions. The relationship between the 
change in variance and the diffusion coefficient is 
given by 

D = - (10) 

2 At 


The diffusion coe f f ici cn r is equal to one half the change 
in variance (An 2 ) divided by the time interval (At). 

By dividing equation 9 by the maximum concen- 
tration at the centroid ( w max ) * taking the log of each 
side and multiplying by 2, the equation becomes 



Le 1 1 i n g 



equation 11 reduces to that of an ellipse. 


X 2 
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I! 

b 2 
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( 12 ) 


(13) 


1 


(14) 
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where a and b are the major and minor semi axis of an 

ellipse fitted to a line of equal concentration about 

the dye patch. Since the edge of the dye patch generally 

formed an irregularly shaped boundary and is characterized 

by relatively steep concentration gradient, the visible 

boundary of the patch was assumed to have a concentration 

of W /2 . With this assumption equations 12 and 13 
max 

reduce to 


0 . 72a e 2 

(IS) 

0.72b 2 
c 

(16) 


By computing the variances in both the longitudinal and 
transverse direction from two flights over the area, the 
diffusion coefficients can be determined from equation 10. 

Since the dye patches seldom form a perfect 
ellipse, the major and minor semi axes of the patch can 
not be measured directly. In the computations, the 
irregular-shaped dye patch is replaced with an equiva- 
lent ellipse that has the same horizontal area and the 
same ratio of transverse to longitudinal variances as 
the real patch. The photographic coordinates of both 
the control points and the outlines of tne dye patches 
are measured with an x-y coordinatograph . The control 
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points are necessary for photographic orientation as 
described in the following section. 

PHOTOGRAPH ORIENTATION 

Photographic orientation is accomplished by 
a non-linear solution to the collinearity condition 
equations (Keller and Tewinkel, 1966). Corrections 
arc generally not required for atmospheric refraction, 
earth curvature, film shrinkage or lens distortion. 

The relationship between photo coordinates and ground 
coordinates is 


x - X 


" X - Y 1 

p o 


p c 

y p - ^o 

= K [RM] 

Y - Y 
P c 

- f 


1-1 



p c _ 


where x and y p are photo coordinates of image point p, 
f is the camera focal length, x Q and y Q are photo co- 
ordinates of the principal point, K is a scale factor, 
the X, Y, and Z subscripted p and c refer to ground 
coordinates of the object and camera station respectively, 
and RM is the rotational matrix. The matrix is defined as 
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"mil ml2 ml3‘ 

[RM] = m21 m22 m23 (18) 

_m31 m32 m33_ 

where: mil = cos(<}>) cos (K ) 

a 

ml2 = cos(W 0 ) sin(K ) + sin(W 0 ) sin(i}>) cos(K ) 

a a 

ml3 = sin(W 0 ) sinfK ) - cos(W 0 ) sin(4>) cos(K ) 

a a 

m21 = -cos(<J>) sin(K 0 ) 

mZ2 = cos(W 0 ) cos(K ) - sin(W 0 ) sin($) sin(K ) 

3 3 

m23 = sin(Wo) cos(K ) + cos(W 0 ) sin(<}>) sin(K ) 

m31 = sin(4>) 

m32 = sin(W 0 ) cos((J>) 

m3 3 = cos (Wo) cos(<j)) 

The throe parameters W 0 , <J> and K are the photographic 
rotations about the X, Y and Z axis, respectively. 

The collinearity equations are obtained by 
dividing the first and second "ows of equation 17 by the 
third row hereby eliminating the scale factor. 

Xp - Xo mil (X p - X c ) + ml2 (Y p - Y c ) + ml3 (Z p - Z £ ) 

-f = m31 (X p - X c ) + m3 2 (Y - Y c ) + m33 (Z - l Q ) 

(19) 

y p - y 0 _ m21 (Xp - X c ) + m2 2 (Y p - i c ) + m23 (Z p - Z Q ) 

Tf m31 (X p - X c ) + m32~(Y^~ : Y c ) + m33 (Z p - Z Q ) 

( 20 ) 
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These equations insure that the camera station, image 
and object lie on a straight line. For each point two 
collinearity equations can be written. As there are 
six unknowns (X - Y . Z - W 0 , $, and K„) a minimum of 
three noncollinear control points are required for 
their solution. However, a least-squares solution 
permits the use of an unlimited number of control points. 

The solution to the equations is obtained based 
on a set of initial approximations which are adjusted 
iteratively until the adjustments become small. The 
collinearity equations are linearized by the Taylor 
series with the expansion terminated at the first de- 
rivative . 

When the initial approximations of B ^ are 
close to the actual parameters values (B) 

f(B) = F (B j ) + E AB i (21) 

Letting 


Y = f (B) - f(B..) 

7 _ 9 f (B) 

1 ~ 9B i 


( 22 ) 


for B = B j 


(23) 
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then 


6 

y = I z. AB. + e 
i= 1 1 1 


which is a linear form of the collinearity equations. 
The least squares solution in matrix notation is 


B = [Z T Z] _1 Z 1 ^ 


(24) 


The initial approximations of the parameters (B^) are 
replaced by 


B j+1 = B. + AB . (25) 

This iterative process is continued until the solution 
converges, that is, until all AB's are less than some 
prespecified amount. In this space resection problem 
no test is made on the linear adjustments but the solu- 
tion is terminated when the angular adjustments are 
less than about two seconds of arc. 

By knowing the photo orientation, it is pos- 
sible to determine the position vectors for any point 
on the photograph. The position vector (X) based on 
the state plane coordinate system axis is related to 
the photographic vector (Xp) by the equation 
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X= [RM] “ 1 Xp (26) 

TIic ground coordinate of any point on the 
photograph can be computed from the unit position vec- 
tor and the camera station coordinates. 


COMPUTATIONS 


The coordinates of the outline of the dye 
patches were measured with the x-y coordinatograph for 
both flights. Photo coordinates about the dye patch 
are converted to state plane coordinates and the cen- 
troid and moments of inertia computed from the fallow- 
ing equations . 


n 

Z-vi 

i=l 


< Y i-l 


- Vi> 


5T = 

Y - 


I 


x 




AA i 



(27) 

(28) 

(29) 

(? 0 ) 
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V ' L V “i ( 3I > 

r xy - D X i Y i “i t32) 

Care must be taken when programming these equations to 
avoid roundoff error. State plane coordinates are typi- 
cally six to seven digit numbers. 

The rotation angle (a) of the axis to obtain 
the principal axis is given by 


Tan (2a) 
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xy 


! x ’ X y 


(3 3) 


The maximum and minimum moments of inertia about the 
principal axis are computed from 


[ = Ix -~ . !y - + \ C+: ly\ + ia 

max 2 \ / xy_ 




(34) 


i . - LLLi . 

mm 7 


[M • ■■ J ' 


/* 


(35) 


The irregular shaped dye patch is replaced 
with an equivalent ellipse that has the same horizontal 
area. In addition, the ratio of I max to I m ^ n are the 
same for the ellipse as for the dye patch. The principal 
moments of inertia for an ellipse are given by 
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* max 


I . 
min 


Trba, 3 
e e 


na 3 
e e 


(36) 


(37) 


where a and b are the major and minor semi axes of an 
equivalent ellipse. If the ratio (HI) of I to I 
remain the same for both the dye patch and the ellipse, 
the ratio of the variances are also the same (equations 
15 and 16) . 


IR 


I 

max 

I . 
mm 




(38) 


After computing the area (A) and the ratio of principal 
moments of inertia (IR) for the dye patch, the semi axes 
of the equivalent ellipse are determined from 

a 2 = - (IR) - 1 / 2 (39) 

e 7T 

b 2 = - (IR)' 1 / 2 (40) 

C „ 

Equations 15 and 16 provide the relationship 
between the semi axes and the variances of the equivalent 
ellipse. The real dye patch is replaced with an ellipse 
that has the same area and the same ratio of transverse 


ak.o-.-1 
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to longitudinal variances as the irregular-shaped, 
original patch. The diffusion coefficients are deter- 
mined from the change in variances by equa' ion 10 while 
the current velocity is determined, from the change in 
position of the centroid of the dye patch between flights. 
The required time interval between photographic flights 
depends on the accuracy of orientation and positioning, 
but generally 30 to 60 minutes is adequate. 

SUMMARY 


This report presented a simplified procedure 
for determining water current velocities and diffusion 
coefficients. Dye drops which form dye patches in the 
receiving water are made from an aircraft. The changes 
in position and size of the patches are recorded from 
two flights over the area. The simplified data proces- 
sing procedure requires only that the ground coordinates 
about the dye patches be determined at the time of each 
flight. With an automatic recording coordinatograph for 
measuring coordinates and a computer for processing the 
data, this technique provides a practical method of de- 
termining circulation patterns and mixing charac.teristi cs 
of large aquatic systems. 
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An effective waste management program for the 
preservation of water quality requires that the fate and 
effect of waste loadings to be predicted before they are 
allowed to be discharged. Mathematical models of the re- 
ceiving waters can be used for making these predictions; 
however, they require as input the mixing and circula- 
tion characteristics of the aquatic system. Remote 
sensing offers a feasible method of measuring both the 
diffusion coefficients and the water currents of a re- 
ceiving waters. This information is necessary to assess 
the environmental impact of waste water discharges and 
for industrial plant siting. 
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